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I. 


INTRODUCTION 


Description  of  the  scattering  of  both  acoustic  and  electro¬ 
magnetic  waves  by  localized  turbulent  structures  (turbules)  in  the 
atmosphere  is  currently  of  great  interest.  The  usual  model  of  an 
individual  turbule  consists  of  a  localized  region  in  which  the 
atmospheric  mass  density  differs  from  that  of  the  surroundings; 
such  a  model  can  be  described  in  terms  of  a  variable  refractive 
index.  (See,  for  example,  Ishimaru  1978) .  The  principal  purpose 
of  the  work  reported  here  was  to  investigate  the  effect  on  acoustic 
scattering  of  including  realistic  flow  velocities  in  a  model 
turbule;  to  the  best  of  the  author's  knowledge,  this  has  not  been 
done  previously.  In  order  to  achieve  this  purpose,  it  was  first 
necessary  to  develop  the  appropriate  acoustical  wave  equations  and 
descriptions  of  the  scattering  from  the  fundamental  fluid 
equations. 

This  report  is  organized  as  follows.  In  section  II,  the 
above-mentioned  theoretical  structure  that  was  developed  is  pre¬ 
sented,  along  with  the  appropriate  scattering  theory.  The  basic 
fluid  equations  were  chosen  to  be  those  for  an  ideal  fluid  with 
negligible  heat  flow.  The  principal  results  quoted  in  this  section 
are  i)  that  the  logarithm  of  the  ratio  of  the  fluid  particle  number 
density  inside  a  quasistatic  turbule  to  that  of  the  surroundings  is 
proportional  to  the  square  of  the  ratio  of  the  turbule  flow  speed 
to  the  acoustic  wavespeed  in  the  surroundings,  in  the  absence  of 
sound  waves,  and  ii)  the  acoustical  wave  equations,  in  terms  of 
small  space  and  time-dependent  variations  in  number  density  and 
fluid  velocity  with  respect  to  their  local  quasistatic  values, 
include  terms  that  are  first  order  in  this  ratio.  Since  this  ratio 
is  generally  quite  small,  the  presence  of  such  first-order  terms  in 
the  wave  equations  presages  the  possibility  that  the  contributions 
of  these  terms  to  the  scattering  might  dominate  those  of  the 
density  variations. 

In  section  III,  the  application  of  the  theory  to  a  realistic 
quasistatic  nonuniformly  spinning  turbule  is  presented.  Two 
approaches  to  calculating  the  acoustical  scattering  are  discussed. 
One  is  the  Born  approximation.  Analysis  is  given  to  show  that  this 
approximation  should  be  reasonably  valid  for  turbule  size  parame¬ 
ters  ka  <  5  and  very  good  for  size  parameters  <  1.  (Here,  k=2ir/ 
(incident  wavelength),  a  =  effective  radius).  Another  is  the 
digitized  Green  function  (DGF)  method,  which  should  be  applicable 
to  a  turbule  of  arbitrary  structure  and  size  (Goedecke  and  O'Brien, 
1988) .  However,  during  the  course  of  the  work  reported,  the  author 
was  unable  to  apply  this  method  to  the  complete  equations;  its 
application  is  thus  far  restricted  to  the  equations  in  which  only 
the  first-order  terms  in  the  above-mentioned  ratio  are  retained. 
The  principal  results  quoted  in  this  section  are  the  Born  approxi¬ 
mation  scattering  efficiency  vs.  scattering  angle  for  several 
values  of  size  parameter,  and  of  total  scattering  efficiency  vs. 
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size  parameter.  The  major  findings  are  that  i)  the  above-mentioned 
first-order  terms  in  the  wave  equations  produce  scattering  cross 
sections  about  three  orders  of  magnitude  larger  than  do  the  terms 
involving  density  changes  alone,  for  realistic  values  of  the  flow 
parameters,  over  a  range  of  size  parameters  from  0.25  to  2.0,  and 
ii)  these  first-order  efficiencies  are  proportional  to  (ka)6,  while 
the  efficiencies  due  to  density  variability  have  the  expected  (ka)4 
dependence,  for  ka  <  2. 

In  section  IV,  a  brief  discussion  is  given  of  the  results  of 
the  work  performed  and  of  needs  for  future  development.  In 
particular,  it  is  suggested  that  even  more  realistic  quasistatic 
turbule  models  should  be  investigated.  These  include  "smoke 
rings,"  i.e.,  toroidal  flow  patterns.  But  this  requires  another 
development  from  first  principles,  because  such  vortices  must 
translate  as  a  whole  at  constant  speed  in  still  air  in  order  to 
exist.  Also,  it  is  suggested  that  the  DGF  method  be  developed 
further  and  eikonal  methods  also  be  developed  for  application  to 
the  problem. 
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II.  THEORY 


A.  Fluid  Equations 

1.  General.  We  consider  a  fluid  with  negligible  viscosity  and 
heat  conduction.  For  simplicity,  we  let  the  fluid  be  a  single¬ 
component  ideal  gas;  extension  to  a  multi-component  fluid  is 
straightforward.  The  fluid  ^equations  involve  the  number  density 
n(?,t),  the  velocity  field  (V(r,t),  the  pressure  p(?,t),  the  tem¬ 
perature  T(? ,t) ,  the  internal  energy  density  U(r,t),  and  the  mass 


the  atoms,  as 

follows: 

dn/dt  +  V-(nV) 

=  0 

(i) 

av/at  +  v-w  = 

— ♦ 

-\^)/mn 

(2) 

au/at  +  v*(vu) 

+  p  V-V  =  0 

(3) 

p  =  nkBT  =  TU, 

r  ■  Y"i  / 

(4) 

Y  =  cp/cv  =  ratio  of  specific  heats 

The  first  two  equations  are  those  of  number  (mass)  continuity 
and  the  Euler  equation  for  an  ideal  fluid.  The  third  is  the 
thermal  energy  flow  equation;  it  results  from  the  more  complete 
energy  continuity  relation 

0t  (^nmV2  +  U)  +  V-[V(U+p)  +  3smnV2V  +  Q]  =  0,  (5) 

when  the  heat  flux  Q  =  0,  in  view  of  eqs.  (1)  and  (2) .  Eq.  (4)  is 
the  ideal  gas  equation  of  state  and  the  expression  of  equipartition 
for  the  model  fluid. 

Combining  (1),  (3),  (4)  yields 

T'1  (0T/0t  +  V'VT)  +  rV-V  =  0  (6) 


Rewriting  (1)  yields 


rn'1  (dn/0t  +  V'Vn)  +  rV*V  =  0 


(7) 


Combining  (6)  and  (7)  yields 


0t  [T  «n(n)  -  £nT]  +  V-V[r«n(n)  -  £nT]  =  0 


(8) 
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Since  we  have  assumed  zero  heat  flux  (adiabatic  conditions) ,  we 
should  have 

p  =  Cn7,  c  =  constant  (9) 

But,  since  p  =  nkBT  also,  we  have 

T  =  Cnr/kB  (10) 

But  then  eq.  (8)  is  satisfied  identically.  Thus,  we  can  eliminate 
T  from  eqs.  (l)-(4).  In  particular,  we  may  write  (2)  as 

dv/dt  +  V-W  =  -  (ykBT/m)  n_1Vn  =  -  (yC/rm)Vnr  (11) 

We  assume  that  the  fluid  contains  a  localized  region  R  cen¬ 
tered  at  the  origin^of  coordinates  in  which  there  is  a  quasistatic 
departure  of  n  and  V  from  their  uniform  constant  background  values 
at  spatial  infinity.  Also,  we  assume  that  there  can  be  a  small- 
amplitude  wavelike  disturbance  in  n  and  V  everywhere  in  the  fluid. 
That  is,  we  assume 

n(r,t)  =  nQ(?) (1  +  e(r,t)),  (12) 

V(?,t)  =  v(r)  +  u(r,t) ,  (13) 

where  |e|  «1  and  |u|  is  "small."  Also,  we  assume  that  as  r-«», 

n0(r)  -*  n*  =  const,  v(r)  -*•  0.  (14) 

The  form  of  eqs.  (7)  and  (11)  suggests  that  we  define 

£n (n/n*)  a  *(r)  +  e(r,t),  (15) 

where  we've  used  (12)  and  kept  terms  to  first  order  in  e(r,t)  only, 
and  ♦(?)  is  defined  by 

*(f)  =  £n(n0(r)/n.)  (16) 

As  we  shall  see  below,  |f(r)|«l  for  realistic  turbules.  So, 
when  we  substitute  (12,  13,  15)  into  the  fluid  equations  (1)  and 
(11),  keep  terms  of  order  (e,u,f,  v2,  cf,uf)  only,  and  require  that 
the  time-dependent  terr..s  and  the  static  terms  in  the  equations  be 
separately  equated,  we  get 


v’Vv  =  -  c2„Vf 

(17) 

v-Vf  +  V-v  =  0 

du/dt  +  v-Vu  +  u-Vv  =  -  c2,V ( A e ) 

(18) 

de/dt  +  V’Vc  +  u*V^r  +  V*u  =  0, 

8 


where 


i  +  r* 

(19) 

=  ykaTa/m 

(20) 

is  the  square  of  the  asymptotic  wavespeed.  Equations  (17)  are  the 
static  (zeroth  order  in  e,u)  equations  that  a  quasistatic  turbule's 
motion  and  density  should  satisfy,  for  adiabatic  conditions. 
Equations  (18)  are  the  acoustic  equations  for  a  wavelike  distur¬ 
bance  in  the  presence  of  the  turbule. 

Now,  as  r-*oo,  f=0,  v=0,  so,  if  Vxv  =  0,  we  have  from 

♦  =  -  v2/2c2„  (21) 

For  cases  of  interest,  v/c0)<<1,  so  |t|  «<1;  ♦  is  second  order 
in  the  small  ratio  v/c„.  If  Vxv  *  0,  which  is  likely,  it  is  still 
clear  from  (17)  that  |f|  is  of  order  v^c2.,  for  if  v  =0  everywhere, 
then,  in  view  of  the  boundary  conditions,  f  a  o  everywhere. 

In  the  acoustical  equations  (18) ,  we  may  consider  time 
dependence  e  iut,  since  these  equations  are  linear  in  (e,u)  and  none 
of  the  coefficients  depend  on  time.  Thus,  we  get 

-iwu  +  Ci  V(Ae )  +  v-Vu  +  U*Vv  =  0  (22) 

-iwe  +  V-u  +  v*Ve  +  u-Viji  =  0 

If  we  are  given  if(r)  and  v(r) ,  localized  in  R,  and  given  a 
plane  wave  incident  from  infinity,  then  in  principle  we  can  solve 
these  equations  for  (e,u)  inside  R,  and  then  predict  how  the 
incident  wave  scatters. 

2.  Wave  equations.  The  first  step  in  obtaining  wave  equations  of 
the  appropriate  form  for  describing  scattering  is  to  decouple  eqs. 
(22)  .  This  is  algebraically  difficult,  but  becomes  reasonably 
simple  if  we  consistently  discard  terms  of  order  v3  and  higher, 
remembering  that  f  is  of  order  v2.  In  the  ensuing  manipulation,  it 
helps  immensely  if  we  put  the  equations  in  dimensionless  form.  To 
do  this,  we  define  or  redefine  several  quantities,  as  follows: 

k  =  o/c„  =  2tt f \ ,  X  =  incident  wavelength 

r  -*  kr  a  r;  V-*  kV(kr)  *  kV  (23) 

v/c.  -♦  v;  u/c.  -*  u 
S  =  Ae 
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With  these  substitutions,  and  keeping  only  terms  up  to  0(v2)  ,  eqs. 
(22)  become 

u  =  -iV<5  +  F 

S  -  -i  [AV*u  +  v'VS  +  u-Vf]  (24) 

where 

F  *  -  i  (v-Vu  +  u'Vv)  (24a) 

It  is  easiest  to  eliminate  S  from  eqs.  (24).  After  some  algebra, 
the  following  equation  results: 

(^+1)5  =  ( V®+1)  F  +  V(U  -  V-F)  a  S(r)  (25) 


where 


Us-  u-Vi|i  -  r^V-u  -  iu*v  +  iv’F  (26) 

Eq.  (25)  is  in  the  "standard  form"  that  allows  a  Green's  function 
treatment  of  the  scattering  problem;  the  source  function  S  goes  to 
zero  rapidly  for  r>a. 

B.  Scattering  Theory 

1.  General.  We  consider  a  plane^wave  of  unit  amplitude  for  the 
relative  density  disturbance  fi(r)  incident  from  infinity  in 
direction  k.  Then 

eln(r)  =  exp(ik-r)  -  uln(?)  =  k  exp(ik-f)  (27) 

This  form  for  uin(?)  follows  from  eqs.  (22)  with  t  s  0,  v  =  0.  Now, 
we  write  down  the  Green's  function  implicit  solution  for  eq.  (25), 
for  outgoing  scattered  waves: 

u(r)  =  uin(r)  -  (47T)'1  /  d3r' (eiR/R)S(r  •)  (28) 


where 


R  a  | r  -  r*  |  (29) 

and  the  integral  on  r'  is  over  all  space. 

The  second  term  on  the  RHS  of  (28)  is  the  scattered  wave 
u,  (?)  .  For  r  -*  «>,  we  get 

u,(r)  (keir/r) f f (k, f )  (30) 

where  the  scattering  amplitude  f(k,f)  is 
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(31) 


f(k,f)  =  -^wkj'^-jd^’  exp(-if-r’)S(r') 

The  form  (30)  of  the  scattered  wave  for  r  -*»  follows  because  u5(?) 
must  also  satisfy  eqs.  (21)  for  $=0,  v=0;  so  we  must  have  for  r-*» 

u,(r)  =  iVe,(r)  =  fc,( r) .  (32) 

This  says  that  the  outgoing  spherical  scattered  wave  must  be  longi¬ 
tudinal  . 

It  is  straightforward  to  show  that  the  standard  definition  of 
the  differential  scattering  cross  section,  equal  to  the  average 
(acoustical)  power  scattered  per  unit  solid  angle  divided  by  the 
incident  power  per  unit  area,  leads  to  the  identification  of  the 
relevant  cross  sections  just  as  in  standard  scattering  theory  in 
quantum  mechanics  or  electromagnetism. 

In  the  scattering  of  electromagnetic  waves  by  a  spherical 
dielectric  particle  of  radius  a,  the  important  quantities  are  the 
size  parameter,  ka  =  2na/X,  and  the  scattering  efficiencies,  which 
are  equal  to  the  scattering  cross  sections  divided  by  the  geometri¬ 
cal  cross  section  ira2.  In  our  treatment  of  the  acoustical  problem, 
an  effective  radius  a  will  be  a  parameter  of  any  model  turbule. 
Thus,  the  important  scattering  efficiencies  will  be  the  following: 

Q(k,f)  =  |  f  (k,f)  |2/ira2  =  differential  scattering  efficiency 

(33) 

Q,(k)  =  JdnQ(k,f)  =  total  scattering  efficiency  (34) 

Q.(k)  =  (4/ka2)  Im(f(k,k))  =  extinction  efficiency  (35) 

Eq.  (35)  is  part  of  the  "optical  theorem,"  which  states  that 
Q„  =  Qa  +  Qa »  Qa  =  absorption  efficiency.  Here,  Qe  =  Qs,  because  we 
have  not  included  viscosity,  so  there  is  no  absorption. 

If  we  substitute  eq.  (25)  for  S  into  eq.  (31)  ,  we  get 
immediately  on  integration  by  parts 

f(k,f)  =  -i(47rk)'1  Jd3r’  (U(r ' ) -V”F (r ' )  )exp(-if-r  •)  .  (36) 

Unfortunately,  we  cannot  simply  evaluate  this  integral,  because  the 
integrand  contains  the  unknown  wave  u(r')  and  its  derivatives.  In 
what  follows,  we  first  derive  the  Born  approximation  ~  t  this 
problem;  then  we  derive  a  digitized  Green’s  function  approach 
(Goedecke  and  O’Brien,  1988) . 

2.  First  Born  Approximation.  This  consists  of  replacing  the 
actual  wave  u(r')  by  the  incident  wave  uin(?')  wherever  the  former 
occurs  in  the  integrand  of  (36) .  In  the  corresponding  problem  in 
electromagnetic  scattering  by  dielectric  particles  of  average 
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radius  a,  refractive  index  m,  this  approximation  is  valid  only  for 
|m-l|ka  «  1.  We  don't  have  a  true  "refractive  index"  in  this 
problem;  but  in  what  follows  we  will  specify  conditions  under  which 
this  approximation  should  be  valid  for  this  acoustic  scattering 
problem. 

We  have  treated  only  cases  for  which  V*v  =  0,  so  what  follows 
is  valid  only  for  this  condition.  After  considerable  algebra,  we 
get  from  (24a,  26,  36) 

f  =  f(1)  +  f(2,n  +  f(2)v  (37) 


where 


f(1>  =  -  (2rrk)‘1  cosSfil^i 


f(2)n  =  +(4rrk)‘1(l-Y  +  l-cos6)I<2)n  (38) 

f(2>v  =  -  (47rk)-1{(l-cos0)  (I(2)v)u 
where  summation  convention  is  used,  and 


I(1)i  ■  Jd3r '  vA  (r ' )  exp(iK-r') 

I(2,n  a  Jd3r'4r(r')exp(iK-r')  (39) 

(I<2,v)ij  3  %/d3r'vi(r')vJ(r')exp(iK-r') 


where 


K  ■  k-f  (40) 

Here,  f(1)  is  the  contribution  that  is  first  order  in  the 
turbule  velocity  v;  f<2)v  is  second  order  in  this;  and  f<2,n  is  the 
only  contribution  that  would  occur  if  the  turbule  density  were 
nonuniform  but  the  flow  v  were  set  equal  to  zero. 

This  is  as  far  as  we  can  go  without  specifying  a  flow  v(r)  and 
a  logarithmic  density  variation  t(?)  for  the  turbule.  We  report 
this  in  the  next  section. 

3.  Digitized  Green  Function  Solution.  This  method  was  originally 
developed  to  allow  modelling  of  electromagnetic  scattering  by 
dielectric  particles  of  arbitrary  shape  and  morphology  (Goedecke 
and  O'Brien  19C8)  .  The  method  would  normally  begin  from  the  causal 
Green's  function  solution  (28)  of  (25).  In  order  to  use  the  DGF 


12 


method,  we  would  have  to  reduce  the  Green's  function  implicit 
solution  first  to  the  form 

Mr)  =  uln>i(r)  +  Jd^'K^r,?*^?')  (41) 

where  is  a  known  function  involving  v,  i|f,  and  the  Green's  func¬ 
tion,  and  their  derivatives.  We  have  not  yet  been  able  to  do  this. 

We  have  been  able  to  express  a  Green's  function  solution  for 
e(r),  valid^only  to  terms  of  first  order  in  v,  in  this  form.  The 
result  for  V*v  =  0  is 

e(r)  =  €ln(r)  +  Jd3r  'K(r,r ' )  e  (? ' )  (42) 


where 


K  =  Ka  +  Kb  (43) 

KA  ( ? ,  r ' )  -  (47T)  1  e^MiR-'-R-2)  (2iVi(r ')  +  8  •  sd  •  jVi  (r) )  (44) 

KB  ( r , r ' )  =  (47t)  _1  eiR(— iR"1  — 3iR”3+3R"3)RiRj(dj'vi(f ' )  )  (45) 


This  result  is  usable,  since  we  know  that  v«l  in  application,  so 
the  first  order  terms  in  v  should  usually  dominate  second  order 
contributions  to  the  scattering. 

In  order  to  proceed  with  the  DGF  method,  we  have  to  choose 
analytic  expressions  for  Vi(r'),  so  that  we  may  evaluate  the  deri¬ 
vatives  d'jVif?'),  dj'dj'Vifr')  at  any  point  r'.  Then  K (?,?')  is  a 
known  function.  Then  we  subdivide  the  region  R  into  cubical  cells 
of  some  side  length  d,  with  cell  number  a  centered  at  r  =  ?a.  Then 
eq.  (45)  becomes  a  matrix  equation, 

^ot  —  ^in.a  ^a/3  ^  0  (46) 

in  obvious  summation  notation,  where 

-  d3K(ra,r*)  .  (47) 

Then  we  invert  this  matrix  equation  numerically,  given  ein  a  =  exp 
( ik'fa) ,  and  obtain  values  for  the  ea,  taking  care  with  the  self¬ 
terms  K^.  Once  the  ea  are  known,  we  can  evaluate  the  scattering 
amplitude  easily,  and  then  obtain  the  cross  sections  equally 
easily.  For  this  formulation,  a  little  algebra  yields  for  the 
scattering  amplitude 

f(k,f)  =  fA  +  f8  (48) 
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where 


fA(k,f)  =  id3(47rk)’1fi  Zaexp  (-if-ra)  (2iVi(rJ  +V2vi(ra))ea 


fB(k,f)  =  id3(4irk)'1f-1fj  Zaexp(-if-ra)  (d^r,,) )  ea 

Experience  with  the  DGF  method  in  electromagnetic  scattering 
has  shown  that,  in  order  to  achieve  accurate  results,  the  cell  side 
length  d  must  be  chosen  <  |  m  | _1 ,  where  m  =  refractive  index  of  scat- 
terer.  (Remember,  here  4  is  in  units  of  (incident  wave  length/ 2  tt)  )  . 
If  the  average  linear  extension  of  a  turbule  is  a  (length  units), 
then  the  number  of  cells  needed  is  >(|m|ka)3.  It's  not  clear  here 
just  what  to  use  for  an  effective  refractive  index;  but,  roughly 
speaking,  the  total  wavespeed  at  any  point  £  in  the  turbule  lies 
between  c„„(l  +  v(r)))  and  c.(l-v(r)).  Thus,  the  effective  refrac¬ 
tive  index  satisfies  l-v(r)  <  m,£f  <  1  -  v(r)  .  Since  v(r)  «  1  in 
applications  here,  then,  the  number  of  cells  needed  for  accuracy 
with  the  DGF  method  must  be  >  (ka)3.  Clearly,  for  large  size 
parameters  ka,  the  inversion  of  (46)  can  be  very  demanding  of 
computer  storage  and  time. 
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III.  APPLICATIONS  AND  RESULTS 


A.  Model  Turbule. 

We  have  considered  a  nonuniformly  spinning  turbule  model.  We 
chose  the  following  function  for  the  turbule  flow, 

v (?)  =  (3xr)  exp(-r2/  (ka)2)  (49) 


where 


a  s  =  &va(ka)'x,  va  a  fla/c*,  (50) 
so  that  n  is  an  angular  velocity  parameter. 

The  Gaussian  is  similar  to  what  has  been  used  recently  for 
density  variations  in  a  turbule  with  v=0  (McBride,  1989) .  The 
length  parameter  a  is  a  measure  of  the  "radius"  of  the  turbule. 

Eqs.  (17)  should  be  satisfied.  However,  we  have  been  unable 
to  find  any  spatially  limited  solutions  of  the  first  of  these  equa¬ 
tions  for  (i»,v)  .  With  the  above  v,  we  used  the  result  that  would 
obtain  if  Vxv  were  equal  to  zero, 

*  =  -VV2  -  -  h  a2(r2-(d*r)2)exp(-2r2/  (ka)2) ,  (51) 

as  a  simple  choice  for  ♦  which  satisfies  continuity  and  has  the 
correct  magnitude,  but  does  not  satisfy  the  first  of  eqs.  (17). 
The  results  should  be  indicative  of  the  magnitude  and  angular 
dependence  of  the  scattering  cross  sections  of  actual  turbules,  as 
long  as  this  model  is  quasistatic,  i.e.,  as  long  as  this  density 
variation  and  velocity  structure  persist  for  a  time  that  is  long 
compared  to  the  period  of  the  incident  wave.  We  assume  that  this 
is  the  case.  We  note  that  the  above-mentioned  density  variations 
alone  that  have  been  used  for  a  (quasistatic)  turbule ^cannot  per¬ 
sist  either,  according  to  eq.  (17),  for  if  v=0,  then  Vi|r=0. 

B.  Born  Approximation 

1.  Expressions  for  Scattering  Amplitude.  We  substitute  the 
velocity  v  and  disturbance  If  of  eqs.  (49  and  51)  into  Eqs.  (39)  . 
After  some  algebra,  the  contributions  to  the  scattering  amplitude 
of  eqs.  (38)  are  found  to  be 

fU)(k,f)  =  -  (a/2ffk)cos8eJklfJfikJi  (52) 
f<2>n(k,f)  =  (a2/4trk)(l-y  +  K2/2)%(«tJ  -  M jJJij  (53) 
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f<2,v(k,f)  =  (a2/4wk)  [  (K2/ 4 )  (tfij-fii&j)  +  krM.Vr.iC.Pj]Jij 


(54) 


where 


JA  s  i%ff3/2(ka) 5  Ki  exp(-(ka)2K2/4)  (55) 

Jij  =  h (ka) 5  [Su  -  3j(ka)2  KiKJ](rr/2)3/z  exp(-(ka)2K2/8)  (56) 

and,  as  in  eq.  (40) , 

K,  *  ^  ,  K2  -  2  ( l-cos0)  (57) 

where  0  is  the  polar  scattering  angle,  such  that  k*f  =  cos0.  Note 

that  the  observation  direction  is  specified  by  the  unit  radial 
vector  f, 

t  =  ex  sin0cos0  +  e2  sin0cos0  +  e3cos0  (58) 


in  terms  of  the  angle  0  and  the  azimuthal  angle  <p,  where  the  are 
the  Cartesian  unit  basis  vectors. 

With  no  loss  of  generality,  we  may  take  the  unit  incident  wave 
vector  to  be  in  the  z-direction, 

k  =  e3  (59) 

The  orientation  of  the  spin  axis  of  the  model  turbule  is  then  arbi¬ 
trary,  so  we  have  the  unit  spin  axis  vector  given  by 

&  =  e3  sinOaCos^,,  +  e2sin0asin0a  +  e3  cos0a  (60) 

In  practice,  it  is  easy  to  see  that  all  differential  cross 
sections  will  depend  on  the  difference  (<p-<pa)  only.  Thus,  we  can 
put  <pa  —  0  during  calculation,  and  then  simply  replace  <p  by  <p-<pa  in 
the  results. 

As  discussed  earlier,  the  (first)  Born  approximation  should  be 
valid  if  the  effective  refractive  index  m  satisfies  |m-l|ka  «1. 
Here,  the  "effective  index"  can  probably  be  taken  as  m  «  1  ±  v„  as 
discussed  in  the  previous  section  on  the  DGF  method.  Thus,  the 
criterion  for  validity  is 

v,ka  «  1  (61) 


16 


In  practice,  |m-l|ka  up  to  unity  is  often  used.  For  realistic 
models,  va  <  0.2,  so  the  Born  approximation  results  should  be 
fairly  trustworthy  for  ka  <  5,  and  very  good  for  ka  <  1,  and  better 
the  smaller  va. 

2.  First  Order  Efficiencies.  If  we  combine  eqs.  (33,  34,  50,  52, 
55,  and  57-60) ,  we  obtain 

Q(1,(kf)  =  (v2a/ 16)  (ka) 6sin20cos20sin20asin2 (0-0a)  (62) 

X  exp[-(ka)2(l-cos0) ] 

This  has  rather  striking  dependence  on  the  observation  direction 
and  on  the  turbule  spin  axis  orientation.  It  goes  to  zero  if  i) 
0a  =  0  or  n ,  ii)  <p  =  <pa  or  <pa  ±  n,  iii)  0  =  0,  n/2,  n.  For  exam¬ 
ple,  suppose  that  0O  =  tt/2,  4>a  =  0.  Then  the  scattering  is  zero  at 
at  0  =  ( 0 , 7t )  and  0  =  (0,^/2,  v) ,  and  is  max  at  <p  -  (n/ 2,  3n/2)  and 
0  =  (rr/4 ,  3tt/4)  .  For  ka  >  1,  the  exponential  becomes  important, 
and  strongly  reduces  scattering  in  the  backward  hemisphere  relative 
to  the  forward. 

This  differential  efficiency  also  has  an  unusual  dependence  on 
ka,  going  as  (ka)6,  quite  different  from  Rayleigh  scattering,  which 
goes  like  (ka)4  with  very  different  dependence  on  scattering  angle. 

It  is  of  great  interest  to  consider  the  efficiency  averaged 
over  random  orientations  of  the  spin  axis.  Since 

(4rr)'1/dna  sin2  0a  sin2(0-$tt)  =  1/3  (63) 

We  get 


<Qa) (k, f ) >  =  (v*a/48)  (ka)6  sin20cos20  exp[-(ka)2(i-cos0)  ]  (64) 

This  is  plotted  vs.  0  for  several  different  values  of  ka  in  Figs. 
1  and  2,  for  va  =  0.1. 

The  total  scattering  efficiencies  are  obtained  by  integrating 
(62  and  64)  over  all  observation  angles.  This  can  be  done  analyti¬ 
cally.  After  some  algebra,  we  get 

<Qfl,s(k)>  =  (7tv2./6)  (ka) 6  C  bnIn((ka)2)  (65) 

where 

bt  =  0.5,  b2  =  -  1.25,  b3  =  1,  b4  =  -  0.25  (66) 

and 

In(x)  s  Jo2<H  £n  e'*?  (67) 
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avg  scattering  efficiency  (xIO 


beta=  1.0  va— 0.1 


scattering  angle  theta 


Figure  1.  Scattering  efficiency  vs.  scattering  angle  for  a  spin¬ 
ning  turbule,  including  only  first  order  contributions 
for  small  size  parameters. 


beta=1.0  vo=0.1 


scattering  angle  theta 


Figure  2.  Scattering  efficiency  vs.  scattering  angle  for  a  spin¬ 
ning  turbule,  including  only  first  order  contributions 
for  large  size  parameters. 
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A  plot  of  Qa),(k)  vs.  ka  is  given  in  Fig.  3,  for  v, 


0.1. 


beta=1.0  va=0.1 


Figure  3.  Average  scattering  efficiency  vs.  size  parameter  for  a 
spinning  turbule,  including  only  first  order  contribu¬ 
tions. 


It  is  important  to  look  at  limits  for  ka«l  and  ka  large.  For 
example,  for  ka  <  0.1,  In((ka)2)  «2n*V  (n+1)  .  Thus,  for  ka  <  o.l, 
Qa)„(k)  is  proportional  to  (ka)6.  However,  for  (ka)  >  4,  ln((ka)2) 
«  nl/jka)2*1*2,  so  the  behavior  of  Qa),(k)  is  given  by  that  of  the 
leading  term  in  the  sum  (65) , 

Qa,s (k)  ka>4  (wvz./12)  (ka)2.  (68) 

This  proportionality  to  (ka)2  is  typical  of  the  Born  approximation, 
and  is  wrong  for  very  large  ka.  Just  as  in  electromagnetic  scat¬ 
tering,  the  total  scattering  efficiency  should  be  independent  of  ka 
for  very  large  ka. 

3.  Second  Order  Efficiencies.  If  we  combine  eqs.  (33,  34,  50,  52, 
53,  54,  56,  and  57-60),  we  get  expressions  for  the  second  order 
scattering  efficiencies.  These  are  much  too  lengthy  and  compli¬ 
cated  to  write  down  in  detail.  The  method  is  given  in  Appendix  A. 
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In  particular,  we  calculated  only  the  averages  over  spin  axis 
orientation.  The  general  form  is 

<Q(2)  (k,  f )  >  =  (v\/512)  x2  e'*f  E*.0  Fn(x)C"  (69) 

where 

x  s  (ka)2/2  ,  C  =  l-cos0  (70) 

and  the  Fn(x)  are  complicated  polynomials  in  x  and  x2.  It  is  these 
functions  that  are  very  lengthy  expressions.  But,  for  x  «  1,  we 
find  simply 

Fo(x)  =  (1-Y)2  +  (46/ 3 ) ( 1-y )  +  86/15 

FJx)  =  2(1+6)  (1-y  +  26/3)  (71) 

F2(x)  =  ( 1+6) 2 

F3(x)  =  F4  (x)  =  0, 

where  we  have  dropped  all  terms  containing  (x,  x2)  .  Here,  6  is 
either  zero  or  unity.  If  we  choose  6=1,  we  get  the  second  order 
result  that  includes  both  f<2,n  and  f(2)v,  eqs.  (53)  and  (54)  .  If  we 
choose  6=0,  we  get  the  second  order  result  for  a  density  variation 
only,  where  the  i|r  function  is,  however,  still  given  by  eq.  (51)  . 

The  expression  (69)  for  <Q2(k,£)>,  using  complete  expressions 
for  the  Fn(x) ,  is  plotted  in  Figs.  4  through  7  against  scattering 
angle  0,  for  several  values  of  ka,  for  v,  =  0.1,  for  y  =  5/3.  Note 
that  its  maxima  are  about  two  orders  of  magnitude  smaller  than 
those  of  the  first  order  contributions,  as  we  would  expect  for 
small  v,  and  for  these  values  of  ka. 

The  total  scattering  efficiency  is  then  given  by  the  integral 
over  solid  angle  of  (69) .  The  result  is 

<Q<2)  (k)  >  =  (ttv\/256)x2  Z‘.0  Fn(x)In(x)  (72) 

where  the  In(x)  are  defined  by  eq.  (67).  This  scattering  effi¬ 
ciency  is  plotted  as  a  function  of  ka  in  Figs.  8  and  9.  We  note 
from  eqs.  (72)  and  (67)  that,  for  small  ka  (small  x)  ,  this  is 
proportional  to  (ka)*,  not  (ka)6  as  for  the  first  order  contri¬ 
bution.  This  means  that,  for  large  enough  v.  and  small  enough  ka, 
the  second  order  contributions  could  dominate  the  first  order.  But 
for  the  realistic  parameter  choices  made  here,  namely  v,  =  0.1  and 
0.25  <  ka  <  2,  the  first  order  is  totally  dominant. 

Note  that  there  is  no  interference  between  first  order  and 
second  order  contributions.  Eqs.  (52)  through  (57)  reveal  that  the 
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first  order  scattering  amplitude  is  purely  imaginary,  while  the 
second  order  ones  are  purely  real.  Thus,  when  we  write 

f  =  f(i>  +  f(2>  =  i|fu>|  +  |  f<2>  |  ,  (73) 

then  we  get  simply 

!  f  j  2  =  |f(1)|2  +  |  f<2>  |2  (74) 

Thus,  the  principal  contributions  of  the  flow  to  the  efficiencies 
go  like  v2,,  while  the  leading  contributions  of  the  density 
variation  go  like  v*a. 
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IV.  DISCUSSION 


The  principal  theoretical  and  calculational  results  of  the 
work  reported  are  those  mentioned  in  the  introduction.  We  rephrase 
them  somewhat  here:  i)  In  order  for  a  quasistatic  density  varia¬ 
tion  to  exist  in  an  isentropic  ideal  fluid  in  the  absence  of  gravi¬ 
tational  and  Coriolis  effects,  a  quasistatic  flow  must  also  exist, 
and  according  to  eqs.  (17)  the  relative  density  variation  is  pro¬ 
portional  to  the  square  of  the  small  ratio  of  flow  speed  to 
asymptotic  acoustical  wavespeed,  ii)  The  fundamental  fluid  equa¬ 
tions  yield  wave  equations  (18)  for  small-amplitude  acoustic 
disturbances  superimposed  on  the  quasistatic  flow  and  density 
variation  of  a  model  turbule.  These  wave  equations  contain  terms 
linear  in  the  flow  speed  ratio,  which  should  then  produce  more 
scattering  than  the  associated  density  variations;  and  iii)  The 
first  Born  approximation  calculations  of  scattering  by  a  non- 
uniformly  spinning  Gaussian  model  turbule  show  that  the  contri¬ 
butions  to  the  scattering  of  the  terms  linear  in  flow  speed  in  the 
wave  equations  are  indeed  about  three  orders  of  magnitude  larger 
than  the  contributions  of  the  terms  involving  the  density  varia¬ 
tions,  for  size  parameters  from  0.25  to  2.0,  for  flow  speed  ratios 
<  0.1.  For  example,  Figs.  3  and  8  show  that  the  first  order  total 
scattering  efficiency,  due  to  the  associated  quasistatic  density 
variations  alone,  is  about  10‘7,  for  a  size  parameter  of  unit,  and 
a  flow  speed  ratio  of  0.1. 

Further  investigation  seems  desirable,  and  is  planned. 
Efforts  probably  should  include  i)  Realistic  modeling  of  other 
quasistatic  turbules  that  include  flow,  such  as  smoke  rings 
(toroidal  vortices) .  This  would  require  modification  of  the 
scattering  equations,  since  such  vortices  must  translate  in  still 
air;  ii)  Extension  of  the  DFG  method;  iii)  Possible  application  of 
eikonal  approximations  similar  to  those  that  have  been  used  with 
great  success  for  large  size  parameters  in  quantum  mechanics  and, 
recently,  in  electromagnetic  scattering  by  dielectric  particles 
(Chen,  1989) . 
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Second  order  Born  approximation 


beta=0.0  va=0.1 


Figure  4.  Scattering  efficiency  vs.  scattering  angle  for  small 

size  parameters,  for  turbule  with  density  variation  only 
(no  flow) . 


Second  order  Born  approximation 


scattering  angle  theto 


Figure  5.  Scattering  efficiency  vs.  scattering  angle  for  large 
size  parameters,  for  turbule  with  density  variation 
only  (no  flow) . 
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Second  order  8orn  approximation 


Figure  8 


Figure  9 


beta=0.0  vo=0.1 


Average  scattering  efficiency  vs.  size  parameter  for 
the  case  of  density  variation  only  (no  flow)  in  the 
turbule. 


Second  order  Born  approximation 
beta®  1.0  va=0.1 


Average  scattering  efficiency  vs.  size  parameter  for 
a  spinning  turbule,  including  only  second  order 
contributions. 
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APPENDIX  A.  ORIENTATION-AVERAGED  CROSS  SECTIONS 


The  total  second  order  scattering  amplitude  of  eqs.  (53)  and 

(54)  , 

f<2)  _  fn<2)  +  fv<2),  may  be  written  as 


f(2)  =  (a2/87rk)  J (x,0  F(x,£,q>,0a)  (Al) 

where 

£  3  1  -  cos0 ,  x  =  (ka) 2/2  (A2) 

J(x,£)  *  rr3/2  x  5/2  e'x{  (A3) 

F(x,s,<p,ej  =  c(x,o  +  s(x,$,(p,e0)  (A4 ) 

C(x,0  =  A(5)  (2-xZ)  (A5) 

A(0  =  %[l-Y+Ul+fi)  3  (A6) 


S(x,£,<p,0a)  =  h  x  A(£)D(£,<p #0a)  +  Bsin20a(l-%xsin20sin2<p) 

(A7) 

D(5/Vf0a)  =  £2cos20a  +  sin20sin20acos2«p-2£sin0sin0acos0acosq> 

(A8) 

We  want  to  calculate  <F2>,  where 

<G(0a/<p)>  s  (47 r)’1  f0 2*  dq>  //  d0a  sin0a  G(0a,  <p  ) 

(A9) 

with  G(0a,<P)  any  functions  of  these  variables. 

We  see  that 

<F2>  =  C2  +  2C<S>  +  <SZ> 

(A10) 

In  order  to  evaluate  the  various  averages,  we  need 
<cos?>  =  <sin<p>  =  0,  <cos2f>  =  <sin2<p>  =  \ 

<sin20a>  =  2/3,  <cos20a>  =  1/3  (All) 

<sin*0a>  =  8/15,  <cos*0a>  =  3/15 

<sin20acos20a>  =  2/15,  <cos*f>  *  <sin4*»  =  3/8 
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From  these,  we  get 

<D>  =  25/3  ,  <D2>  =  452/5 

<Dsin26a>  =  85/15  -  252/15  (A12) 

<Dsin20asin2<p>  =  25/15 
Then ,  we  get 

<S>  =  (26/3)  [l-X(25-52)/4]  (A13 ) 

<S2>  =  x2A252/5  +  86/15  -  (46/15)  (25“52)X  (A14) 

+  6(25-52)2  x2/20  +  6Ax(85/15  -  252/ 15-5  (25“52)  x/15) 

Then 

2C<S>  =  6(2-x5) (4A/3) [l-X(25-52)/4]  (A15) 

C2  =  A2(2-x5)2  (A16) 

If  we  combine  these  according  to  eq.  (A10)  ,  we  may  rearrange  terms 
so  that  the  result  looks  like 

<F2>  =  S*n«0  F„(X)  5“  (A17) 


Using  (A2,A1)  and  a  =  (v,)/(ka)  ,  x  *  (ka)2/2,  yields  the 

expression  (69)  of  the  text. 
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APPENDIX  B.  SUMMARY  OF  RELEVANT  FORMULAS 


A  model  spinning  turbule  was  chosen  to  have  the  following 
quasistatic  flow  velocity  v(?)  and  molecular  number  density 

n0(f): 

v(?)  =  (fixf)exp  (-r2/a2)  (Bl) 

n0(f)  =  -njj2  f  (2c2)  (B2) 

This  represents  a  non-uniformly  spinning  turbule,  with  angular 
velocity  parameter  n  about  an  axis  n  through  the  origin  of  coor¬ 
dinates.  The  quantities  (n<,,c«)  are  the  background  (number  density, 
wavespeed)  far  from  the  turbule;  the  parameter  a  is  an  "effective 

radius"  of  the  turbule;  r  is  the  position  vector.  The  expression 

for  n0(f)  follows  approximately  from  the  assumption  that  turbules 
are  formed  adiabatically . 

If  an  acoustic  plane  wave  with  propagation  vector  ic  ,  where 
k  *  2n/\  ~  w/ c„,  o  *  angular  frequency,  is  incident  upon  this  model 
turbule,  then  the  scattering  amplitude  f(  £  ,t)  in  first  Born 


approximation  is  given  by 

=  f<x)  +  +  42)  (B3) 

where  f(1)  is  linear  in  fi,  and  f(2)  is  quadratic  in  n.  The  differen¬ 
tial  scattering  cross  section  then  turns  out  to  be  given  by 

o(£,f)  =  |f<l>  I2  +  |42)  +  f'2)|2  ■  0{1)  +  0(2)  (B4) 

The  expressions  for  f(1)  and  f(2)  are  as  follows: 

f (1)  (£,  f)  =  -  (a2/2nk)  cos0eJJCi  AkJi  (B5) 

=  (a2/4rcic)  (1  -  Y  +  K2/2)'/2(bij  -  fljO,)  (B6) 
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(B7) 


fl2)  (£,f)  =  (a2/4itic){(/C2/4) 

+  M.<W  rJi 


where 


a  =  Q/o  (B8) 

JA  =  i  %  tt3/2  (ka)5  Ki  exp(-(ka)2K2/4)  (B9) 

Jtj  =  ^(ka)5  [Sij  -  ^(ka)2KiKj]  ( w/ 2 ) 3/2  exp (- (ka) 2K2/8  (BIO) 

and 


Kk  S  kt  -  ft,  K2  =  2(1  -  cos  6)  ,  (Bll) 

when  0  is  the  polar  scattering  angle,  the  angle  between  k  and  the 

observation  direction  $  -  ^sinBcos^  +  £2sin0sin<J>  +  e3cos0.  These 

equations  (B5)  through  (Bll)  are  the  same  as  eqs.  (50)  and  (52) 
through  (57)  of  the  text.  Here  <p  is  the  azimuthal  scattering 
angle. 

The  differential  cross  section  is  given  by  inserting  (B5) 
through  (B7)  into  (B4)  .  It  depends  on  the  azimuthal  angle  <p  as 
well  as  on  the  polar  angle  0.  For  example,  the  scattering 
efficiency  from  a1*3  =  |  f a)  | 2  is 

Q{1)  (k,  f )  =  a(1)/ica2  ,B:m 

=  (Qa/4c.) 2  (ka) 6  [sin0cos0sin0Qsin  (<j>  -  <j)Q)  ] 2  1  ' 

where  (0fl/0n)  are  the  (polar,  azimuthal)  angles  of  the  spin  vector 

n. 


There  is  interest  in  the  average  of  these  differential  cross 
sections  over  random  orientations  of  the  spin  axis.  These  were 
obtained  in  two  ways,  analytic  and  numerical.  The  analytic  method 
is  outlined  in  the  text  and  in  Appendix  A.  The  numerical  method 
simply  numerically  averages  over  uniformly  spaced  values  of  cos0n 
and  0ft.  The  results  of  this  latter  method  agreed  with  those  of  the 
former,  which  are  displayed  as  Figs.  1  through  9  of  the  text.  The 
analytic  relation  is 

(Qn)(k,f))-  (2/3)  (Qa/icJ2  (ka)6  (sin0cos0)2  (B13) 
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If  Ei  is  the  incident  exitance  (watt.m-2)  of  the  sound  field  and 
1(0)  is  the  intensity  (watt.sr'1)  of  the  remote  scattered  field  at 
the  off-axis  angle  0,  then 

I (0)  =  xa2(s?(1)  (k,t))E1  (B14) 

The  cross  section  a<2)  is  much  smaller  than  a(1),  for  realistic 
choices  of  the  parameter  fia/c.  «1. 

The  total  cross  sections,  defined  by 

o (i)  (k,  f )  =  fdQoU)  (k,t)  ,i  =  1,2  (B15) 

or  rather  the  scattering  efficiencies  Q(i)  =  aa)/Tia2,  were  also 
obtained  both  analytically,  as  described  in  the  text  and  in 
appendix  A,  and  numerically,  by  numerical  integration  over  (0,0). 
Here,  d  n=  sin  0  d  0d  <p  =  solid  angle  element.  The  results  agreed. 
The  analytic  total  scattering  efficiency  for  a(1>  is 

<0a)(ic)>=  U/45)  (Qa/cJ 2  (ka) 6  (B16) 


Some  authors,  especially  in  radar  scattering,  relate  the  differen¬ 
tial  scattering  cross  section  to  the  ideal  isotropic  scatterer. 
The  scattered  field  from  the  ideal  isotropic  scatterer  is  the  same 
for  all  scattering  angles  0  with  magnitude  of  the  total  cross 
section  times  the  incident  exitance  divided  by  4 n.  The  actual 
field  from  a  non-isotropic  scatterer  is  then  the  magnitude  of  the 
previous  sentence  times  an  appropriate  function  of  0  called  the 
phase  function.  This  information  is  provided  to  alert  the  reader 
that  there  is  a  quantity  of  4rr  applied  in  the  radar  scattering 
concept  that  is  not  needed  in  the  development  of  this  report. 
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TEXT  MODE  REPRESENTATIONS 


GRAPHICS  MODE  REPRESENTATIONS 


n 

n 

n 

n 

ft 

A 

A 


Q  OMEGA  bold 

H  OMEGA  underline 

A  OMEGA  hat 

fi  OMEGA  vec 
fl  OMEGA  bar 


A  A 

Q  Q 

r  tuvWERx tuvWER 


stack  {‘#n> 


q  stack  {A#n> 

^  stack  {A#n> 

^  stack  {/Vfl} 

q  func  stack  {A#n> 


stack  {func  A#  func  n} 
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